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(57) Abstract: An apparatus 10,70 and 
method is provided that includes a spatial 
array of unsteady pressure sensors 15-18 
placed at predetermined axial locations Xi 
-Xn disposed axiaJly along a pipe 14. For 
measuring at least one parameter of a solid 
particle/fluid mixture 12 flowing in the pipe 
14. The pressure senors 15-18 provide 
acoustic pressure signals Pi(t)-pN(t) to a 
signal processing unit 30 which determines 
the speed of sound a^ix of the particle/fluid 
mixture 12 in the pipe 14 using acoustic 
spatial array signal processing techniques. 
The primary parameters to be measured 
include fluid/particle concentration, 
fluid/particle mixture volumetric flow, and 
particle size. Frequency based sound speed 
is determined utilizing a dispersion model 
to determine the parameters of interest. 
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APPARATUS AND METHOD FOR MEASURING PARAMETERS OF A MIXTURE HAVING SOLID 
PARTICLES SUSPENDED IN A FLUID FLOWING IN A PIPE 

Cr OSS-Reference To Related Patent Applications 

5 This application claims the benefit of U.S. Provisional Application No. 60/35 1,232 

(ridra Docket No. CC-0410), filed January 23, 2002; U.S. Provisional Application No. 
60/359,785 (Cidra Docket No. CC-0403), filed February 26,2002; U.S. Provisional 
Application No. 60/375,847 (Cidra Docket No. CC-0468), filed April 24, 2002; U.S. 
Provisional Application No. 60/425,436 (Cidra Docket No. CC-0538), filed November 12, 
10 2002; and U.S. Provisional Application No. 60/426,724 (Cidra Docket No. CC-0554), filed 
November 1 5, 2002, all of which are incorporated herein by reference in their entirety. 

Technical Field 

This invention relates to an apparatus for measuring the flow passing Avithin a pipe, 
15 and more particularly to an apparatus and method for measuring the speed of sound and/or 
vortical disturbances propagating in the flow, having particles suspended within a 
continuous fluid, to determine parameters, such as particle/fluid ratio, particle size and 
volumetric flow rate of the flow in pipes using acoustic and/or dynamic pressures. 

# 

20 Background Art 

This invention provides a method to measure parameters of a fluid/particle mixture 
in a pipe that can be used in many applications, such as in chemical, pharmaceutical, 
petroleum and power generation industries. Li particular, the invention provides a method 
to measure pulverized coal and air mixtures used in pulverized fiiel delivery systems in 

25 place in a large percentage of coal fired boilers used in the power generation, industry. 

. Currently, well over 50% of the electricity in the US is generated with coal. While 
coal is considered a cost effective, abimdant resource in the US, the use of coal has been 
restricted due in large part to environmental concerns. To mitigate this impact, the US 
Department of Energy and the Power Generation industry have large programs designed to 

30 develop technology to reduce the environment efiTects of burning coal. These Clean Coal 
hiitiatives include technology designed to develop improvements in the combustion process 
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to improve efificiency while reducing pollutants such as unbumed carbon, ash, and nitrous 
oxide (NOx). 

The ability to measure the flow rate and composition of the air /coal mixture within 
the coal pipes is an important aspect of any system or strategy designed to optimize the 
5 performance of the PF delivery system. The industry recognizes this and therefore has been 
developing a wide variety of technologies to perform this measurement. These include 
probe based and sampling devices, as well as real time meters based on a wide variety of 
technologies including electrostatic charges, microwaves, and ultrasonic. 

1 0 Summary of the Invention 

Objects of the present invention include providing a system for measuring the speed 
of sound propagating through a particle/fluid mixtm^e in pipes in coal fired boiler systems 
and related processes, for example, to determine particular parameters of the mixture. 
According to the present invention, an apparatus for measuring at least one 

1 5 parameter of a particle/fluid mixture in a pipe includes a spatial array of at least two 

pressure sensors, disposed at different axial locations along the pipe. Each of the pressure 
sensors measiures an unsteady pressure within the pipe at a corresponding axial location. 
Each of said sensors provides a pressure signal indicative of the unsteady pressure within 
the pipe at said axial location of a corresponding one of said sensors. A signal processor, 

20 responsive to said pressure signals, provides a signal indicative of the at least one parameter 
of the mixture in the pipe. 

According to the present invention, a method for measuring at least one parameter of 
a particle/fluid mixture in a pipe includes measuring unsteady pressures within the pipe at at 
least two predetermined axial measxurement locations along the pipe to provide a pressure 

25 signal indicative of the unsteady pressure within the pipe at each of the at least two 

predetermined axial measurement locations. Further the method includes calculating the at 
least one parameter of the particle/fluid mixture in the pipe using the imsteady pressure 
measured at the axial measurement locations. 

The foregoing and other objects, features and advantages of the present invention 

30 will become more apparent in light of the following detailed description of exemplary 
embodiments thereof. 
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Brief Description of the Drawings 

Fig. 1 is a block diagram of a flow meter for measuring the speed of sound of the 
fluid/particle mixture flowing with a pipe, in accordance with the present invention. 
5 Fig. 2 is a schematic diagram of a pulverized fuel (PF)/air mixture parameter 

measurement system within a coal fired boiler system, in accordance with the present 
invention. 

Fig. 3 is a magnified photograph showing particle size of coal typical of the system 
shown in Fig. 2. 

10 Fig. 4 is a plot of the speed of sound of a mixture versus the frequency in air/coal 

mass flow ratio, in accordance with the present invention. 

Fig. 5 is a plot of actual data and a model of the speed of sound as a function of 
frequency for air/coal mixtures, in accordance with the present invention. 

Fig. 6 is a plot showing the standard deviation of sound speed versus frequency for 
15 various arrays of a PF/air mixture parameter measurement system, in accordance with the 
present invention. 

Fig. 7 is a plot of sound speed as a function of frequency for air/coal nfiixtures. with 
fixed particle size (50 mm) and varying air-to-fuel mass Ratio in accordance with the 
present invention. 

20 Fig. 8 is a plot of sound speed as a function of frequency for air/coal mixtures with 

varying particle size where the air-to-fuel mass ratio is equal to 1.8 in accordance with the 
present invention. 

Fig. 9 is a plot of soimd speed as function of air/coal ratio in accordance with the 
present invention. 

25 Fig. 10 is a flow diagram of an optimization procedure employed to determine air- 

to-fiiel ratio and particle size from analytical model and experimentally determined 
dispersive speed of sound data in accordance with the present invention. . 

Fig. 1 1 is a plot of the results of the optimization procedure of Fig. 10 applied to 
data recorded from an array of sensors listening to flow in a six inch circular duct, 50 \im 

30 particle size, 100 fl/sec air flow rate with an air-to-fiiel ratio of 1.8. 
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Fig. 12 is a plot of the results of the optimization procedure of Fig. 10 applied to a 
series of data sets with varying air-to-fuel ratio. 

Fig. 13 is a ko) plot of data processed from an array of pressure sensors use to 
measure the speed of sound of a coal/air mixture flowing in a pipe, in accordance with the 
5 present invention. 

Fig. 14 is a block diagram of a flow meter for measuring the vortical field of the 
fluid/particle mixture flowing with a pipe, in accordance vnih the present invention. 

Fig. 15 is a cross-sectional view of a pipe showing a turbulent pipe flow velocity 

profile. 

1 0 Fig. 1 6 is a side elevational view of another embodiment of a flow meter for 

measuring the vortical disturbances in a pipe, in accordance with the present invention. 

Fig. 17 is a plot of the pressure signals measured by a pair of pressure sensors of the 
flow meter of Fig. 16. 

Fig. 18 is a plot of the cross-correlation of the pressure signals plotted in Fig. 17. 
15 Fig. 19 is a ko) plot of data processed from a flow meter embodying the present 

invention that illustrates slope of the convective ridge, and a plot of the optimaztion 
function of flie convective ridge, in accordance with the present invention. 

Fig. 20 is a side elevational view of a plurality of pressure sensors, havmg PVDF, 
clamped to the outer surface of the pipe, in accordance with the present invention. 
20 Fig. 21 is a partial perspective view of one of the pressure sensors of Fig. 20. 

Best Mode for Carrying Out the Invention 

Referring to Figs. 1 and 14, a flow meter 10,70 embodying the present invention is 
provided that measures a mmiber of parameters/characteristics of a mixture 12 of solid 

25 particles suspended within a continuous fluid flowing within a pipe or conduit 14, wherein a 
fluid is defined as a liquid and/or a gas. The flow meter may be configured and 
programmed to measure the speed of sound propagating through the mixture or measure the 
vortical disturbances propagating through the mixture. In some instances, the flow meter 10 
may be configured to measure both the speed of soimd and the vortical disturbances. 

30 Depending on the configuration or embodiment, the flow meter can measure at least one of 
the following parameters of the mixture flow 12: the fluid/particle concentration 
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(volumetric phase fraction), the voliimetric flow rate, the size of the solid particles, the mass 
flow of the mixture and the velocity of the mixture. To determine any one of these 
parameters, the flow meter 10,70 measures the unsteady pressures created by the speed of 
sound (SOS) and/or the vortical disturbances propagating through the mixture flowing in 
5 the pipe 14, 'vhich will be described in greater detail hereinafter. 

The solid particles of the mixture 12 may be of any size, shape and material. For 
example, the particles may be small in size as in the form of a powder, in a granular form, 
or greater in size. The flow meter 10,70 can be used in any application that carries solid 
particles suspended in a fluid through a pipe, such as in chemical, pharmaceutical, 

10 petroleum and power generation applications. For example, the present invention is well 
suited to measure the parameters (e.g. air/coal ratio, particle size) for power generation 
systems that use pulverized coal to fire the furnace a steam boiler system. 

As one example, the present invention will be discussed in the context of a 
Pulverized Fuel (PF) delivery system for power generation, but one will appreciate that the 

15 flow meter can be applied to any number of other applications, as discussed hereinbefore. 
A representative PF delivery system 1 is shown in a coal fired boiler system 2 in Fig. 2. 
The coal is pulverized in a mill 3 and entrained in air produced by many means, such as a 
fan 4 to transport the PF/air mixture via pipes 12 for delivery to the furnace 6. Typical 
fumaces can have >50 coal pipes, each 12-20 inch in diameter. Typically, a large utility 

20 boiler >300 Mw, can have 4-1 1 pulverizing mills feeding the furnace. The ability of the PF 
delivery system to deliver the proper amount of fuel and air to the furnace through these 
multiple coal pipes, both collectively and individually, has a strong influence on the 
performance and emissions from the coal fired boiler. 

As is known, non-uniformities in the PF delivery system 1 can result in variation of 

25 the fiiel to air ratios, causing hot spots, regions of high NOx generation, and unbumed fuel. 
The connection between performance of a PF fiiel delivery system 1 and a boiler system 2 is 
well recognized. The flow meter 10 embodying the present invention is capable of 
measuring the fiiel to air ratio and particle size of the pulverized coal provided to the 
furnace to thereby provide feedback to the operator to provide more efQcient combustion of 

30 the coal. 
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As described hereinbefore, the flow meter 10,70 of the present invention may be 
configured and programmed to measiire and process the detected unsteady pressures Pi(t) - 
PN(t) created by acoustic waves and/or vortical disturbances propagating through the 
mixture to determine parameters of the mixture flow 12. One such flow meter 10 is shown 

5 in Fig. 1 that measwes the speed of sound (SOS) of one-dimensional sound waves 

propagating through the fluid/particle mixture to determine the composition tiie mixture, 
namely the liquid/particle ratio of the mixture. The flow meter is also capable of 
determining the average size of the particles, velocity of the mixture, and the volumetric 
flow rate of the mixture. It is known that sound propagates through various mediums at 

10 various speeds in such fields as SONAR and RADAR fields. The speed of sound of a 

mixture within a pipe 14 may be determined using a number of known techniques, such as 
those set forth in U.S. Patent Application Serial No. 09/344,094, entitled *Tluid Parameter 
Measurement in Pipes Using Acoustic Pressures", filed June 25, 1999, and U.S. Patent 
Application Serial No. 10/007,749, entitled "Fluid Parameter Measurement in Pipes Using 

15 Acoustic Pressures", filed November 7, 2001, each of which are incorporated herein by 
reference. The present invention utilizes at least one flow meter 10 to determine various 
parameters of the liquid/particle mixture, wherein one of the parameters is the speed at 
which sound travels within the mixture pipe system as will be more fully described herein 
below. 

20 In accordance with the present invention, the speed of sound propagating through 

the mixture 12 is measxired by passively listening to the flow with an array of unsteady 
pressure sensors to determine the speed at which one-dimensional compression waves 
propagate through a liquid/particle mixture contained within the pipe 14. 

As shown m Fig. 1, the flow meter 10 has an array of at least three acoustic pressure 

25 sensors 15,16,17, located at three locations xi,X2,X3 axially along the pipe 14. One will 

appreciate that the sensor array may include more than three pressure sensors as depicted by 
pressure sensor 18 at location xn. The pressure generated by the acoustic waves may be 
measured through holes in the pipe 14 ported to external pressure sensors 15 - 18 or by 
o&er techniques discussed hereinafter. The pressure sensors 15-18 provide pressure time- 

30 varying signals Pi(t),P2(t),P3(t),PN(t) on lines 20,21,22,23 to a signal processing unit 30 to 
known Fast Fourier Transform (FFT) logics 26,27,28, 29, respectively. The FFT logics 26 - 
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29 calculate the Fourier transform of the time-based input signals Pi(t) - PnO) and provide 
complex frequency domain (or frequency based) signals Pi(ca),P2(cD),P3(co),PN(co) on lines 
32,33,34,35 indicative of the frequency content of the input signals. Instead of FFT's, any 
other technique for obtaining the frequency domain characteristics of the signals Pi(t) - 
5 PnCO, may be u^ed. For example, the cross-spectral density and the power spectral density 
may be used to form a frequency domain transfer functions (or frequency response or ratios) 
discussed hereinafter. 

The frequency signals Pi(cd) - PhK®) are fed to a^ix-Mx Calculation Logic 38 which 
provides a signal to line 40 indicative of the speed of sound of the mixture amix (discussed 
10 more hereinafter). The amix signal is provided to map (or equation) logic 42, which converts 
amix to a percent composition of the PF/air mixture and provides a %Comp signal to line 44 • 
indicative thereof (as discussed hereinafter). Also, if the Mach number Mx is not negligible 
and is desired, the calculation logic 40 may also provide a signal Mx to line 46 indicative of 
the Mach number Mx. 

15 More specifically, for planar one-dimensional acoustic waves in a homogenous 

mixture, it is known that the acoustic pressure field P(x,t) at a location x along a pipe, where 
the wavelength X of the acoustic waves to be measured is long compared to the diameter d 
of the pipe 12 (i.e., X/d »1), may be expressed as a supeii>osition of a right traveling wave 
and a left traveling wave, as follows: 



where A,B are the frequency-based complex amplitudes of the right and left traveling 
waves, respectively, x is the pressxire measurement location along a pipe, co is frequency (in 
rad/sec, where Q)=27cf), and kr,ki are wave nimibers for the right and left travelling waves, 
respectively, which are defined as: 



where amix is the speed of sound of the mixture in the pipe, a> is frequency (in rad/sec), and 
Mx is the axial Mach number of the flow of the mixture within the pipe, where: 



20 




Eq. 1 



25 




Eq.2 



mix 



Eq. 3 where Vmix is the axial 



^mix 
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velocity of the mixture. For non-homogenous mixtures, the axial Mach nmnber represents 
the average velocity of the mixture and the low frequency acoustic field description remains 
substantially unaltered. 

The data from the array of sensors may be processed in any domain, including the 

5 frequency/spatial domain, the temporal/spatial domain, the temporal/wave-mmiber domain 
or the wave-number/frequency (k-o) domain. As such, any known array processing 
technique in any of these or other related domains may be used if desired. 

Also, some or all of the functions within the signal processing unit 30 may be 
implemented in software (using a microprocessor or computer) and/or firmware, or may be 

10 implemented using analog and/or digital hardware, having sufiBcient memory, interfaces, 
and capacity to perform the fimctions described herein. 

Acoustic pressmre sensors 15-18 sense acoustic pressure signals that, as measured, 
are lower firequency (and longer wavelength) signals than those used for ultrasonic flow 
meters of the prior art, and thus the current invention is more tolerant to inhoraogeneities in 

15 the flow, such as roping and other time and space domain inhomogeneities within the flow, 
even where entrenchment or coal "roping" is unlikely such as following a bend. The term 
"roping" is a term known to those skilled in this art which represents a form of severe 
spatial and temporal mal-distribution induced in mixture flows of widely different 
component densities. It is a condition where a large portion of the coal flow is in a band 

20 mnning along one side of pipe 14. 

In addition, the present invention incorporates the compliance of the pipe 14 to 
detemiine the effective speed of sound of the pipe/PF/air mixture system. The acoustic 
pressure signals Pi(t) - PN(t) are generated within the PF/air mixture of the pipe 14 by a 
variety of non-discrete sources such as remote machinery, mills, fans 4 (Fig. 2), valves, 

25 elbows, as well as the PF/air mixture flow itself It is this last source, the PF/air mixture 12 
flowing within the pipe 14, which is a generic source of acoustic noise that assures a 
minimmn level of acoustics for any PF/air mixture piping systems for which the present 
invention takes unique advantage. The flow generated acoustics increase with mean flow 
velocity and the overall noise levels (acoustic pressure levels) are a function of the 

30 generating mechanism and the damping mechanism. As such, no extemal discrete noise 

source is required within the present invention and thiis may operate using passive listening. 

-8- 



BNSOOCID: <WO ^030e275aAl.L> 



wo 03/062759 



PCTAJS03/01925 



While the flow meter 10 passively listens to the mixture flow 12. the present invention 
contemplates adding an acoustic source to inject a desire acoustic wave into the flow to be 
measured, such as by compressing, vibrating and/or tapping the pipe, to name a few 
examples. 

For certain typ-^s of pressure sensors, e.g., pipe strain sensors, accelerometers, 
velocity sensors or displacement sensors, discussed hereinafter, it may be desirable for the 
pipe 14 to exhibit a certain amount of pipe compliance. 

Alternatively, to minimize any error effects (and the need for the corresponding 
calibration) caused by pipe compliance, the axial test section 50 of the pipe 14 along where 
the sensors 15 - 18 are located may be made as rigid as possible. To achieve the desired 
rigidity, the thickness of the wall of the test section 50 may be made to have a 
predetemnned thickness, or the test section 50 may be made of a very rigid material, e.g.. 
steel, titanium, Kevlar®, ceramic, or other material with a high modulus. 

It is within the scope of the present that the pressure sensor spacing may be known 
or arbitrary and that as few as two sensors are required if certain information is known 
about the acoustic properties of the PF/air mixture piping system. The pressure sensors are 
spaced sufficiently such that the entire length of the array (aperature) is at least a sigmficant 
fraction of the measured wavelength of the acoustic waves being measured. As will be 
described in greater detail, the acoustic wavelength to be measured is a function of at least 
the size and mass of the particles, and the viscosity of the fluid. The greater the size and 
mass of the particles and/or the less viscous the fluid. .the greater the spacing of the sensors 
is needed. Conversely, the smaller the size and mass of the particles and/or the more 
viscous tiie fluid, the shorter the spacing of the sensors is needed. 

As discussed, the flow meter 10 measures the speed of sound of one-dimensional 
sound waves propagating through the fluid/particle mixture to determine the composition of 
the mixtiue. Specifically, the speed of sound propagating through dilute solid/air mixtures 
can be directiy related to the mass fraction particles of the flow. A typical PF fiiel dehvery 
systeml mayoperatewithanairtocoalmassratioof 1.5to2.5. Typically. PF delivery 
systems operate with an air-to-coal.mass ratio of 1 .5 to 2.5 with coal density of 1200 to 
) 1400 kg/m' compared to 1.2 kg/m' for air at standard atmospheric conditions. Hius. 
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meeting the desired mass ratio results in a very dilute mixture of coal on a volumetric basis, 
on the order of one part in 1000 by volume. 

Assuming that the particles of coal are small enough and the acoustic frequencies 
and the frequencies of perturbations associated with the acoustics are low enough for the 
solid particles to exhibit negligible slip ( both steady and unsteady), the sound speed can be 
assumed to be non-dispersive (that is constant with frequency) and the volimietric phase 
fraction of the mixture could be determined through the Wood equation: 



Pmix =S^/Pl 



1 ^ 



15 Including the effect of the compliance introduced by the conduit 12 (in this case a 

circular pipe of modulus E, radius R and wall thickness t) 



1 I ^ 2R 

+ a where cr s - 



Pmix^ measured Pmix^mix 

20 Utilizing the relations above, the speed at which sound travels within the piping 

system of a representative coal / air mixtures is shown in Fig. 4 as a function of air /coal 
mass ratio. For this example, the pure air was assumed to have a density of 1.2 kg/m^3 and 
a sound speed of 365.9 m/s and the coal was assiuned to have a density of 1400 kg/m^3 and 
a sound speed of 2439 m/s. As shown, the effect of increasing coal fraction, i.e. decreasing 

25 air/coal ratio is to decrease the sound speed. Physically, adding coal particles effectively 
mass loads the mixture, while not appreciably changing the compressibility of the air. Over 
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the parameter range of interest, the relation between mixture sound speed and air / coal ratio 
is well behaved and monatomic. 

While the calibration curves based on predictions from first principles are 
encouraging, using empirical data mapping from sound speed to air / coal ratio may result in 
5 improved accuracy of the p/<=^sent invention to measure the air/coal fractions of the mixture. 

However, it has been discovered that the physical properties of pulverized coal/air 
mixtures are generally such that there will be velocity slip at all but very low frequencies 
(on the order of <l-2 Hz for nominally 50 \im coal particles in air). 

Fig. 5 shows the measured speed of sound as a function of frequency for an actual 
10 coal/air mixture 12. The sound speed was measured utihzing passive listening techniques 
of the present invention as described herein. The frequency dependence of the sound speed 
was determined by applymg a Capon array-processing algorithm at multiple narrow 
frequency ranges between 50-300 Hz thereby determining a frequency specific acoustic 
propagation velocity. In this particular example, the data was obtained wherein the coal/air 
15 mixture was flowing at nominally 100 ft/sec with an air-to-coal mass ratio equal to 1.8. The 
coal particles were nominally 50 \xm in size, representative of pulverized coal typically iised 
in power generation and other industrial applications. A magnified view of the coal 
particles that were used for this test is shown in Fig. 3. 

Further shown in Fig. 5, the soimd speed increases with increasing frequency and 
20 asymptotes toward a constant value. The sound speed asymptote at higher frequency is 

essentially the sound speed of air only with no influence of the suspended particles. Also, it 
is apparent fliat the sound speed of the coal/air mixture has not reached the quasi-steady 
limit at the lowest frequency for which sound speed was measured. The sound speed is 
continuing to decrease at the lower frequency limit. An important discovery of the present 
25 invention is that the speed at which sound propagates through dilute particles suspended in a 
continuous fluid is said to be dispersive. As defined herein, flie speed at which acoustic 
waves propagate through dispersive mixtures varies wifli frequency. 

Measuring the sound speed of a mixture 12 at progressively lower and lower 
frequencies becomes mherently less accurate as the total length of the array of pressure 
30 sensors 15-18 (Axaperaturc), which define the aperature of the array, becomes small 

compared to the wavelength of the acoustics. In general, the aperture should be at least a 
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significant fraction of a wavelengtii of the sound speed of interest. In a particular 
embodiment sound speed data was recorded with an array of four sensors, spaced at 12 
inches, for a total aperture of three feet. At 50 Hz, a 1000 ft/sec sound wave has a 
wavelength of 20 ft. Thus, the aperture of this particular array (approx. 36 uiches) spanned 
5 only 3/20ths of a wavelength, and the array's ability to accurately resolve soimd speeds 
below this was clearly impaired. It is an important aspect of the present invention that the 
ability to resolve soxmd speed at low frequencies is directly related to aperture of the array. 
Consequently longer arrays are used to resolve sound speeds at lower frequencies. As 
shown in Fig. 6, the standard deviation associated with detemiining the speed of sound in 
10 air is shown as a function of frequency for three arrays of varying aperture, namely 1.5 ft, 3 
ft and 10 ft. 

Given the practical constraints in accurately measuring soimd speeds at ultra-low 
frequencies, the data suggests that utilizing a quasi-steady model to interpret the 
relationship between sound speed, measured at frequencies above those at which the quasi- 
1 5 steady model is apphcable, and the air-to-fuel ratio would be problematic, and may, in fact, 
be impractical. Thus, the key to imderstanding and interpreting the composition of coal/air 
mixtures through sound speed measurements lies in the dispersive characteristics of the 
coal/air mixture. 

In accordance with the present invention the dispersive nature of the system utilizes 
20 a first principles model of the interaction between the air and particles. This model is 
viewed as being representative of a class of models that seek to account for dispersive 
effects. Other models could be used to account for dispersive effects without altering the 
intent of this disclosure (for example, see the paper titled "Viscous Attenuation of Acoustic 
Waves in Suspensions" by R.L. Gibson, Jr. and M.N. Toksoz), which is incorporated herein 
25 by reference. The model allows for slip between the local velocity of the continuous fluid 
phase and that of the particles. The drag force on the particles by the continuous fluid is 
modeled by a force proportional to the difference between the local fluid velocity and that 
of the fluid particles and is balanced by inertial force: 



30 



dt 
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where K = proportionality constant, Uf = fluid velocity, Up = particle velocity, pp = particle 
density and Vp = particle volume. 

The effect of the force on the continuous fluid phase by the fluid particles is modeled as a 
5 force term in the axial momentui^ equation. The axial momentum equation for a control 
volume of area A and length Ax is given by: 



where P = pressure at locations x and Ax, (t>p = volume fraction of the particles, pf == fluid 
10 density. 

The particle drag force is given by: 



15 

where Cd = drag coefiicient, Ap = frontal area of particle and pf = fluid density. 

Using Stokes law for drag on a sphere at low Reynold's number gives the drag coefficient 
as: 




20 



24 
Re 



24 fi 



where D, 



particle diameter and p = fluid viscosity. 



25 Solving for K in this model yields: 
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Using the above relations and 1-dimensional acoustic modeling techniques, the following 
relation can be derived for the dispersive behavior of an idealized fluid particle mixture. 



In the above relation, the fluid SOS, density (p) and viscosity (</>) are those of the pure phase 
fluid, Vp is the volume of individual particles and <|>p is the volumetric phase fiaction of the 

particles in the mixture. 
10 Two parameters of primary interest in pulverized coal measurements are particle 

size and air-to-fuel mass ratio. To this end, it is of interest to examine the dispersive 
characteristics of the mixture as a function of these two variables. Figs. 7 and 8 show the 
dispersive behavior for coal/air mixtures with parameters typical of those used in pulverized 
coal deliver systems. 

1 5 In particular Fig. 7 shows the predicted behavior for nominaUy 50 nm size coal in 

air for a range of air-to-fuel ratios. As shown, the effect of air-to-fiiel ratio is well defined 
in the low frequency limit. However, the effect of the air-to-fiiel ratio becomes 
indistinguishable at higher frequencies, approaching the sound speed of the pure air at high 
frequencies (above -100 Hz). 

20 Similarly, Fig. 8 shows the predicted behavior for a coal/air mixture with an air-to- 

fiiel ratio of 1 .8 with varying particle size. This figure illustrates that particle size has no 
influence on either the low frequency limit (quasi-steady ) sound speed, or on the high 
frequency limit of the sound speed. However, particle size does have a pronounced effect 
in the transition region. 

25 Figs. 7 and 8 iUustrate an unportant aspect of the present invention. Namely, that 

the dispersive properties of dilute mixtures of particles suspended in a continuous fluid can 
be broadly classified into three frequency regimes: low frequency range, high frequency 
range and a transitional frequency range. Although the effect of particle size and air-to-fiiel 
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ratio are inter-related, the predominant effect of air-to-fuel ratio is to determine the low 
frequency limit of the sound speed to be measured and the predominate effect of particle 
size is to determine the frequency range of the transitional regions. As particle size 
increases, the frequency at which the dispersive properties appear decreases. For typical 
pulverized coal applications, this transitional region begins at fairly low frequencies, ~ 2Hz 

for 50 nm size particles. 

In the low frequency regime, the particles exhibit negligible slip with the fluid. The 
frequency range for which the no-slip, quasi-steady approximation is valid is a fimction of a 
variety of parameters including particle size, continuous phase viscosity, particle shape and 
particle density. 

The quasi-steady sound speed is given by the low frequency limit of the above 
relation, where AFR is air/ftiel ratio: 



i 



1+ 



<PpPp 
Pf 



i 



AFR 



Note that particle size does not affect the low frequency limit of the sound speed. Referring 
to Fig. 9, the sound speed was measured using an embodiment of the present invention 
having eight sensors at 20.5 inch spacing, averaged from 2(M0 Hz, for a range of air-to- 
coal mass ratios. TTie sound speed predicted for the coal/ak mixtures using the quasi-steady 
model are also presented. As shown, although the general trend is captured, i.e. sound 
speed decreases with increased coal loading, the enror is significant, rendering a first 
principle interpretation, based on a quasi-steady model inadequate. 

In the high frequency limit, the dispersion relation predicts the sound speed with 
asymptote towards the sound speed of die pure fluid. 

Interestingly, the high frequency Umit is independent of both particle size and air-to-fuel 
ratio. 
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Given the difficulties measuring sufficientty low frequencies to apply the quasi- 
steady model and recognizing that the high frequency sound speed contains no direct 
information on either particle size or air-to-fiiel ratio, it becomes apparent that the 
dispersive characteristics of tiie coal/air mixture should be utilized to determine particle size 
5 and air-to-fuel ratio based on speed of sound measurements. 

As described hereinbefore, the flow meter 10 of the present invention includes the 
ability to accurately determine the average particle size of the coal in the PF/air mixture 
within the pipe 14 and the air to fuel ratio. Provided there is no appreciable slip between 
the au- and the solid coal particle, the propagation of one dimensional sound wave through 
10 multiphase mixtures is influenced by the effective mass and the effective compressibility of 
the mixture. For an air transport system, the degree to which the no-slip assumption applies 
is a strong function of particle size and frequency. In the limit of small particles and low 
frequency, the no-slip assumption is valid. As the size of the particles increases and the 
frequency of the sound waves increase, the non-slip assumption becomes increasing less 
15 vaUd. For a given average coal particle size, the increase in slip with frequency causes 
dispersion, or, in other words, the sound speed of the mixture to change with frequency. 
With appropriate calibration the dispersive characteristic of a mixture will provide a 
measurement of the average particle size, as well as, the air to friel ratio (particle/fluid ratio) 
of the mixture. 

20 Using the model described above, which yields the equation shown below, and 

experimentally determined sound speed as fimction of frequency, the present invention 
includes an optimization procedure to simultaneously determine particles size and AFR in 
coal / air mixtures: 



25 




Referring to Fig. 10 there is shown an optimization procedure in accordance with the 
present invention in which the fi«e parameters of an analytical model are optimized ti 
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^ini^ize an e„cr Junction. For «Ius«a.ion purposes, 4e «ror function uUliz«. is tte sun. 
of «,c differences of .he sound speeds beWeen an analytical model and toe experunenuny 
detennined sound speed as a function of iirequency: 



5 err- 2-r^^^-^^'"°***^» "^^-^ 
/-/to* 

Tie resulU of 4= optimization procedure applied to data record«l ftom an array of sensors 
listening to flow in a six inch circular duC 50 jun particle size, 100 fVsec air flow rate w,ti> 
anair-tCHfuelratioof 1.8isshowninFig. 11. The measured and opflmized-model- 
,0 predicted sound speed is shown. As shown, the model capBrres the tmnsitional ftequeney 
range well and provides a good estimate of tiie air-to-fuel ratio. 

The results of tire optimization procedure applied to a series of data sets with 
varying air-to-fuel ratio is shown in Fig. 12. Note for this optimization the particle size was 
held constant over the range of data sets. 
,5 ,n addition to measuring the ftaid to particle ratio of the mixtin. 12 and particle Size 

Within a pipe 14 using the measured speed ofsound. the flow meter .0 ^^^"^ 
abiUty to measure of volumetric flow rate of the mixti»e by comparing the drfference of the 
speed of one dimensional sound waves propagating with and agams. the mean flow. 

This method of determining the volumetric flow rate of the particle/fluid maxture 12 
20 within pipe 14 relies on the interaction of the mean flow witi, flte acoustic pressure fleld. 
Tlre interaction results in sound waves propagating with ti.e mean flow traveling at «.e 
3peed ofsound (it the partiCe/liquid mixtirre were no, flowing) plus the convection velocrty 
and. conversely, sound waves traveltag again« the mean flow propagating at ti.. speed of 
sound minus the convection velocity. That is. 



25 



Where a^ = velocity of a right traveling acoustic wave relative to a stationary observer (i.e 
30 thepipe 14).a. = velocity of alefttravelingacousticwaveapparenttoastationary 
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observer, amu - fluid speed of sound (if the fluid were not flowing) and u = the mean flow 
velocity (assumed to be flowing from left to right in this instance). Combining these two 
equations yields an equation for tiie mean velocity, 

•"-^ 

Therefore, by measuring the propagation velocity of acoustic waves in both directions 
relative to the stationary pipe as described hereinbefore, the mean flow velocity can be 
calculated by multiplying the mean flow velocity by the cross-sectional area of the pipe 14. 

The practicality of \ising this method to determine the mean flow is predicated on 
the ability to resolve the sound speed in both directions with sufBcient accuracy to 
determine the volumetric flow. For typical hquid measurements, flow velocities are 
typically at - 10 ft/sec and sound speeds of "-4000 ft/sec. Thus axial mach numbers arc on 
the order of 10/4000 of 0.0025. For a +/- 10% accuracy in flow rate (+/- 1 ft/sec), the sound 
speed of the upstream and downstream propagating waves would need to be resolved to +/- 
0.5/4000 or 1 part in 8,000. 

However, for PF/air mixture flows, axial flow velocities are nominally around 70 
ft/sec with no flow sound speeds of -700 ft/sec. This results in mach numbers of -0.1, 
approximately 2 orders of magnitude greater than typical Uquid flows. For pulverized fuel 
flows, to resolve the flow rate to 10% accuracy (or +/- 7 ft/sec), one would have to resolve 
the sound speed to +/- 3.5 ft/sec, or 3.5/700 or 1 part in 200. 

For the sound speed measurement, the flow meter 10 utilizes similar processing 
algorithms as those employed herein before. The temporal and spatial frequency content of 
sound propagating within the process piping 14 is related through a dispersion relationship. 

k 

fi) = 

The wave number is k, which is defined as k=2n/X, a> is the temporal frequency in rad/sec, 
and amix is the speed at which sound propagates within the process piping. For this cases 
where soimd propagates in both directions, the acoustic power is located along two acoustic 
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ridges, one for the sound traveling with the flow at a speed of a„ux + V^x and one for the 
sound traveling against the flow at a speed of amix - Vmix- 

Fig. 13 shows a k-co plot generated for acoustic sound field of a coal/air mixture 
flowing through a pipe. Two acoustic ridges are clearly evident. Each of the slopes of the 
5 two depicted acoustic ridges respectively defines the speed of sound traveling with and 
against the mean flow. 

The sonar flow meter 10 of Fig. 1 is configured and programmed to measure and 
utilize the speed of sound propagating through a particle/fluid mixture 12 flowing in a pipe 
14 to determine volumetric flow rate. Referring to Fig. 14, a flow meter 70 embodying the 
1 0 present invention includes the ability to measure volumetric flow rate of the mixture by 
measuring the unsteady pressures generated by vortical disturbance 88 propagating in the 
mixture. The flow meter 70 uses one or both of the foUowing techniques to determine the 
convection velocity of the vortical disturbances within the fluid/particle mixture 12 by: 

1) Cross correlating unsteady pressure variations using an array of unsteady pressure 
15 sensors. 

2) Characterizing the convective ridge of the vortical disturbances using an array of 
unsteady pressure sensors. 

20 The overwhelming majority of industrial process flows involve turbulent flow. 

Turbulent fluctuations within the process flow govern many of the flow properties of 
practical interest including the pressure drop, heat transfer, and mixing. For engineering 
applications, considering only the time-averaged properties of turbulent flows is often 
sufficient for design purposes. For sonar flow metering technology, understanding the 

25 time-averaged velocity profile in turbulent flow provides a means to interpret the 

relationship between speed at which coho-ent structures convect and the volumetricaUy 
averaged flow rate within a pipe. 

For turbulent flows, the time-averaged axial velocity varies with radial position, 
from zero at the wall to a maximum at the centerline of the pipe. The flow near the wall is 

30 characterized by steep velocity gradients and transitions to relatively uniform core flow near 
the center of the pipe. Fig. 15 shows a representative schematic of a velocity profile and 
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coherent vortical flow structures 88 present in fiilly developed turbulent pipe flow 12. The 
vortical structures 88 are superimposed ovor time averaged velocity profile within the pipe 
14 and contain temporally and spatially random fluctuations with magnitudes typically less 
than 10% percent of the mean flow velocity. 
5 From a volumetric flow measurement percnective, the volumetrically averaged flow 

velocity is of interest. The volumetrically averaged flow velocity, defined as V=Q/A, is a 
tiseful, but arbitrarily defmed property of the flow. Here, A is the cross sectional area of the 
pipe and Q is the volumetric flow rate. In fact, given the velocity profile within the pipe, 
little flow is actually moving at this speed. 
10 Turbulent pipes flows are highly complex flows. Predicting the details of any 

turbulent flow is one of nature's great-unsolved problems. However, much is known 
regarding the statistical properties of the flow. For instance, turbulent pipe flows contain 
self-generating, coherent vortical structures often termed "turbulent eddies". The 
maximum length scale of these eddies is set by the diameter of the pipe. These structures 
1 5 remain coherent for several pipe diameters downstream, eventually breaking down into 
progressively smaller eddies until the energy is dissipated by viscous effects. 

Experimental investigations have established that eddies generated within mrbulent 
boundary layers convect at roughly 80% of maximum flow velocity. For pipe flows, this 
implies that turbulent eddies will convect at approximately the volumetrically averaged 
20 flow velocity within the pipe. The precise relationship between the convection speed of 
turbulent eddies and the flow rate for each class of meters can be cahTjrated empirically as 
described below. 

The flow meter 70 of Fig. 14 determines the convection velocity of the vortical 
disturbances within the fluid/particle mixfaire by cross correlating unsteady pressure 
25 variations using an array of unsteady pressure sensors, similar to that shown in U.S. Patent 
Application Serial No. 10/007,736. filed November 8, 2001, entitled "Flow Rate 
Measurement Using Unsteady Pressures", which is incorporated herein by reference. 

Referring to Fig. 14, the flow meter 70 includes a sensing section 72 along a pipe 12 
and a signal processing unit 74. The pipe (or conduit) 14 has two measurement regions 
30 76,78 located a distance AX apart along the pipe 14. At the first measurement region 76 are 
two unsteady (or dynamic or ac) pressure sensors 80,82, located a distance X, apart, capable 
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of measuring the unsteady pressure in the pipe 14, and at the second measurement region 
78, are two other unsteady pressure sensors 84,86, located a distance X2 apart, capable of 
measuring the unsteady pressure in the pipe 14. Each pair of pressure sensors 80,82 and 
84,86 act as spatial filters to remove certain acoustic signals from the unsteady pressure 
signals, and the distances Xi,X2 are determined by the desired filtering characteristic for 
each spatial filter, ^ discussed hereinafter. 

The flow meter 70 of the present invention measures velocities associated with 
unsteady flow fields and/or pressure disturbances represented by 88 associated therewith 
relating to turbulent eddies (or vortical flow fields), inhomogeneities in the flow (such as 
bubbles, slugs, and the like), or any other properties of the flow, fluid, or pressure, having 
time varying or stochastic properties that are manifested at least in part in the form of 
unsteady pressures. The vortical flow fields are generated within the fluid of the pipe 14 by 
a variety of non-discrete sources such as remote machinery, pumps, valves, elbows, as well 
as the fluid flow itself. It is this last source, the fluid flowing within the pipe, that is a 
generic source of vortical flow fields primarily caused by the shear forces between the fluid 
and the wall of the pipe that assures a minimum level of disturbances for any fluid piping 
systems for which the present invention takes unique advantage. The flow generated 
vortical flow fields generally increase with mean flow velocity and do not occur at any 
predeterminable jfrequency. As such, no external discrete vortex generating source is 
required within the present invention and thus may operate using passive detection. It is 
wifliin the scope of the present that the pressure sensor spacing may be known or arbitrary 
and that as few as two sensors are required if certain information is known about the 
acoustic properties of the system as will be more fiiUy described herein below. 

The vortical flow fields 88 are, in general, comprised of pressure disturbances 
having a wide variation in length scales and which have a variety of coherence length scales 
such as that described in the reference "Sound and Sources of Sound", A. P.Dowling et al, 
Halsted Press, 1983, which is incorporated by reference to the extend of understanding the 
invention. Certain of these vortical flow fields 88 convect at or near/or related to the mean 
velocity of at least one of the elements within a mixture flowing in a pipe. The vortical 
pressure disturbances 15 that contain information regarding convection velocity have 
temporal and spatial length scales as well as coherence length scales that differ from other 
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disturbances in the flow. The present invention utilizes these properties to preferentially 
select disturbances of a desired axial length scale and coherence length scale as will be 
more folly described hereinafter. For illustrative purposes, the terms vortical flow field and 
vortical pressure field will be used to describe the above-described group of unsteady 
5 pressure fields having temporal and spatial length and coherence scales described herein. 

The pressures Pi,P2,P3,P4 may be measured through holes in the pipe 14 ported to 
external pressure sensors or by other techniques discussed hereinafter/ The pressure sensors 
80,82,84,86 provide time-based pressure signals Pi(t)J>2(t),P3(t), P4(t) on lines 90 - 93, 
respectively, to signal processing unit 74 which provides a convection velocity signal Uc(t) 
10 on a line 96 which is related to an average flow rate Uf(t) of the fluid flowing in the pipe 14. 
Also, some or all of the functions within the signal processmg unit 74 may be 
implemented in software (using a microprocessor or computer) and/or firmware, or may be 
implemented using analog and/or digital hardware, having sufficient memory, interfaces, 
and capacity to perform the fimctions described herein. 
15 In particular, in tiie processing unit 74, the pressure signal Pi(t) on the line 90 is 

provided to a positive input of a sunmaer 100 and the pressure signal PaCt) on the line 91 is 
provided to a negative input of the summer 100. The output of the summer 100 is provided 
to line 104 indicative of the difference between the two pressure signals Pi JPa (e.g.. Pi- 
P2-Pas0. 

20 The pressure sensors 80,82 together with the summer 100 create a spatial filter 76. 

The line 104 is fed to bandpass filter 108, which passes a predetermined passband of 
frequencies and attenxiates firequencies outside the passband. In accordance with the present 
invention, the passband of the filter 108 is set to filter out (or attenuate) the dc portion and 
the high frequency portion of the input signals and to pass the frequencies therebetween. 

25 Other passbands may be used in other embodunents, if desired. Passband filter 108 

provides a filtered signal Pasfl on a line 1 12 to Cross-Correlation Logic 116, described 
hereinafter. 

The pressure signal P3(t) on the line 92 is provided to a positive input of a summer 
102 and the pressure signal P4(t) on the line 93 is provided to a negative input of the 
30 summer 102, The pressure sensors 83,84 together with the sxmimer 102 create a spatial 
filter 78. The output of the sununer 102 is provided on a line 106 indicative of the 
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difference between the two pressure signals P3,P4 (e.g., P3-P4= Pas2). The Une 106 is fed to a 
bandpass filter 1 10, similar to the bandpass filter 108 discussed hereinbefore, which passes 
frequencies within the passband and attenuates frequencies outside the passband. The filter 
110 provides a filtered signal Pasf2 on a line 1 14 to the Cross-Correlation Logic 1 16. The 
5 signs on the summers 1 00, 1 02 may be swapped if desired, provided the signs of both 

summers are swapped together. In addition, the pressure signals P,,P2j*3,P4 may be scaled 
prior to presentation to the summers 100,102. 

The Cross-Correlation Logic 116 calculates a known time domain cross-correlation 
between the signals P^n and T^a on the lines 112,1 14, respectively, and provides an output 
10 signal on a line 11 S indicative of the time delay t it takes for an vortical flow field 88 (or 
vortex, stochastic, or vortical structure, field, disturbance or perturbation within the flow) to 
propagate fix)m one sensing region 76 to the other sensing region 78. Such vortical flow 
disturbances, as is known, are coherent dynamic conditions that can occur in the flow which 
substantially decay (by a predetermined amount) over a predetermined distance (or 
1 5 coherence length) and convect (or flow) at or near the average velocity of the fluid flow. As 
described above, the vortical flow field 88 also has a stochastic or vortical pressure 
disturbance associated with it. In general, the vortical flow disturbances 88 are distributed 
throughout the flow, particularly in high shear regions, such as boundary layers (e.g., along 
the inner wall of the pipe 14) and are shown herein as discrete vortical flow fields 88. 
20 Because the vortical flow fields (and the associated pressure disturbance) convect at or near 
the mean flow velocity, the propagation time delay t is related to the velocity of the flow by 
the distance AX between the measurement regions 76,78, as discussed hereinafter. 

Although pressure disturbances associated with vortical flow fields 88 occur 
naturally in most flow conditions, an optional circumfereaitial groove (not showm) may be 
25 used in the inner diameter of the pipe 14 to help generate unsteady flow fields in the form of 
vertices into the flow. However, the groove is not required for the present invention to 
operate, due to vortex generation which naturally occurs along the pipe inner wall, as 
discussed hereinbefore. Instead of a single circumferential groove a plurality of axially 
spaced circumferential grooves may be used. The dimensions and geometry of the groove(s) 
30 70 may be set based on the expected flow conditions and other factors. Other techniques 
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may be used as vortex generators if desired including those Aat may protrude within the 
inner diameter of pipe 14. 

Referring to Fig. 14, a spacing signal AX on a line 120 indicative of the distance AX 
between the sensing regions 76,78 is divided by the time delay signal x on the line 1 18 by a 
5 divider 122 which provides an output signal on the line 96 i->dicative of the convection 
velocity Uc(t) of the fluid flowing in the pipe 14, which is related to (or proportional to or 
approximately equal to) the average (or mean) flow velocity Ui(t) of the fluid, as defined 
below: 

Uc(t) = AX/xocU<(t) Eq. 1 

10 The convection velocity Uc(t) may then be calibrated to more precisely determine 

the mean velocity Uf(t) if desired. The result of such calibration may require multiplying the 
value of the convection velocity Uc(t) by a calibration constant (gain) and/or adding a 
calibration offset to obtain the mean flow velocity Ui(t) with the desired accuracy. Other 
calibration may be used if desired. For some applications, such calibration may not be 
1 5 required to meet the desired accuracy. The velocities Uf(t),Uc(t) may be converted to 
volumetric flow rate by multiplying the velocity by the cross-sectional area of the pipe. 

Referring to Figs. 16 - 18, as is known, cross-correlation may be used to determine 
the time delay t between two signals yi(t),y2(t) separated by a known distance AX, that are 
indicative of quantities 80 that convect with the flow (e.g., density perturbations, 
20 concentration perturbations, ten^jerature perturbations, vortical pressure disturbances, and 
other quantities), hi Fig. 16, the signal y2(t) lags behind the signal yi(t) by 0. 15 seconds. If 
a time domain cross-correlation is takrai between the two signals yi(t),y2(t), the result is 
shown m Fig. 17 as a curve 124. The highest peak 126 of the curve 124 shows the best fit 
for the time lag x between the two signals yi(t),y2(t) is at 0.15 seconds which matches the 
25 reference time delay shown in Fig. 1 7. 

Referring to Fig. 14, as discussed hereinbefore, since pressure disturbances 
associated withm the vortical flow field 88 convect (or flow) at or near the average velocity 
of the fluid flowing in the pipe 14, the vortical pressure disturbances observed at the 
downstream location are substantially a time lagged version of the vortical pressure 
30 disturbances observed at the iq)stream location. However, the total vortical pressure 

perturbations or disturbances in a pipe may be e:q)ressed as being comprised of vortical 
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pressure disturbances (Pvorticai), acoustic pressure disturbances (Pacoustic) and other types of 
pressure disturbances (Po*er) as shown below expressed in terms of axial position along the 
pipe at any point in time: 

P(x,t) = Pvo«ical (X,t) + Pacoustic (x,t) + Pother (X,t) Eq. 2 

As a result, the unsteady pressure disturbances Pvonica. can be masked by the acoustic 
pressure disturbances Pacoustic and the other types of pressure disturbances Pother- In 
particular, the presence of the acoustic pressure disturbances that propagate both upstream 
and downstream at the speed of sound in the fluid (sonic velocity), can prohibit the direct 
measurement of velocity from cross-correlation of direct vortical pressure measurements. 

The present invention uses temporal and spatial filtering to precondition the pressure 
signals to effectively filter out the acoustic pressure disturbances PacousUc and other long 
wavelength (compared to the sensor spacing) pressure disturbances in the pipe 14 at the two 
sensing regions 76,78 and retain a substantial portion of the vortical pressure disturbances 
Pvomci associated with the vortical flow field 88 and any other short wavelength (compared 
to the sensor spacing) low frequency pressure disturbances Pother- In accordance with the 
present invention, if the low frequency pressure disturbances P<«ber are small, they will not 
substantially impair the measurement accuracy of Pvonicai. 

The Pvortcai dominated signals from the two regions 76.78 are then cross-correlated 
to determine the time delay x between the two sensing locations 76,78. More specifically, 
at the sensing regiori 72, the difiference between the two pressure sensors 80,82 creates a 
spatial filter 76 that effectively filters out (or attenuates) acoustic disturbances for which the 
wavelength X of the acoustic waves propagating along the fluid is long (e.g.. ten-to-one) 
compared to the spacing X, between the sensors. Other wavelength to sensor spacing ratios 
may be used to characterize the filtering, provided the wavelength to sensor spacing ratio is 
sufficient to satisfy the two-to-one spatial aliasing Nyquist criteria. 

Thus, if the pressure sensors P, J>2 have an axial spacing Xi, and assuming that the 
spatial filter 76 will attenuate acoustic wavelengths longer than about 10 times the sensor 
spacing X,. the smallest acoustic wavelength haun that is attenuated would be: 

X™„=10CX,) Eq.3 
One dimensional acoustic disturbances are also governed by the following knovra 
inverse wavelength-frequency relation: 
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X = a/f or f=a/X Eq.4 
where a is the speed of sound of the fluid, f is the frequency of the acoustic 
disturbance, and X is the wavelength of the acoustic disturbance. 

Using Eq. 4, such a spatial filter would filter out frequencies below about: 

5 fmax == a/^min Eq. 5 

The above discussion on the spatial filter 76 also applies to the second spatial filter 
78 comprising the other pair of pressure signals P3,P4, axially spaced a distance X2 apart, 
which provides ttie differenced vortical pressure signal Pas2. 

The second technique of determining the convection velocity of the vortical 
1 0 disturbances within the fluid/ particle mixture is by characterisizing the convective ridge of 
tiie vortal disturbances using an array of unsteady pressure sensors, as will be described. 

The sonar flow metering methodology uses the convection velocity of coherent 
structure with turbulent pipe flows to determine the volumetric flow rate. The convection 
velocity of diese eddies 88 is determined by applying sonar arraying processing techniques 
1 5 to determine the speed at which the eddies convect past an axial array of unsteady pressure 
measurements distributed along the pipe 14. 

The sonar-based algorithms determine the speed of the eddies by characterizing both 
the temporal and spatially fi-equency characteristics of the flow field. For a train of coherent 
eddies convecting past a fixed array of sensors, the temporal and spatial firequency content 
20 of pressure fluctuations are related through the following relationship: 

k 



Here k is the wave number, defined as k=2n/X and has units of 1/length, C0 is the temporal 
firequency in rad/sec, and Uconvcct is the convection velocity. Thus, the shorter the 
wavelength (larger k) is, the higher tfie temporal frequency. 

25 In sonar array processing, the spatial / temporal frequency content of time stationary 

sound fields are often displayed using 'Tc-o plots". K-(0 plots are essentially three- 
dimensional power spectra in which the power of a sound field is decomposed into bins 
correspondmg to specific spatial wave numbers and temporal frequencies. On a k-o plot, 
the power associated with a pressure field convecting with the flow is distributed in regions 

30 which satisfies the dispersion relationship developed above. This region is termed "the 
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convective ridge" (Beranek, 1992) and the slope of this ridge on a k-w plot indicates the 
convective velocity of the pressure field. This suggests that the convective velocity of 
turbulent eddies, and hence flow rate within a pipe, can be determined by constructing a k-co 
plot from the output of a phased array of sensor and identifying the slope of the convective 
ridge. 

Fig. 19 shows an example of a k-eo plot generated from a phased array of pressure 
sensors. The power contours show a well-defined convective ridge. A parametric 
optimization method was used to determine the "best" line representpg the slope of the 
convective ridge 100. For this case, a slope of 14.2 ft/sec was determined. The 
intermediate result of tfie optimization procedure is displayed in the insert, showing that 
optimized value is a unique and well-dejBned optima. 

The k-w plot shown in Fig. 19 illustrates the fundamental principle behind sonar 
based flow measure, namely that axial arrays of pressure sensors can be used in conjunction 
with sonar processing techniques to determine the speed at which naturally occurring 
turbulent eddies convect within a pipe. 

The pressure sensors 15-18 described herein may be any type of pressure sensor, 
capable of measuring the unsteady (or ac or dynamic ) pressures within a pipe 14, such as 
piezoelectric, optical, capacitive, resistive (e.g., Wheatstone bridge), accelerometers (or 
geophones), velocity measuring devices, displacement measuring devices, etc. If optical 
pressure sensors are used, the sensors 15-18 may be Bragg, grating based pressure sensors, 
such as that described in US Patent Application, Serial No. 08/925,598, entitled " High 
Sensitivity Fiber Optic Pressure Sensor For Use In Harsh Environments", filed Sept. 8, 
1997, now U.S. Patent 6,01 6,702. Alternatively, the sensors 14 may be electrical or optical 
strain gages attached to or embedded in the outer or inner wall of the pipe which measure 
pipe wall strain, including microphones, hydrophones, or any other sensor capable of 
measuring the unsteady pressures within the pipe 14. In an embodiment of the present 
invention that utilizes fiber optics as the pressure sensors 14 they may be connected 
individually or may be multiplexed along one or more optical fibers using wavelength 
division multiplexing (WDM), time division multiplexing (TDM), or any other optical 
multiplexing techniques. 



-27- 



wo 03/062759 



PCT/US03/01925 



For any of the embodiments described herein, the pressure sensors, including 
electrical strain gages, optical fibers and/or gratings among others as described herein, may 
be attached to the pipe by adhesive, glue, epoxy, tape or other suitable attachment means to 
ensure suitable contact between the sensor and the pipe 14. The sensors may alternatively 

5 be removable or pemianently attached via known mechanical techninues such as 

mechanical fastener, spring loaded, clamped, clam shell arrangement, strapping or other 
equivalents. Alternatively, the strain gages, including optical fibers and/or gratings, may be 
embedded in a composite pipe. If desired, for certain applications, the gratings may be 
detached from (or strain or acoustically isolated from) the pipe 14 if desired. 

10 It is also within the scope of the present invention that any other strain sensing 

technique may be used to measure the variations in strain in the pipe, such as highly 
sensitive piezoelectric, electronic or electric, strain gages attached to or embedded in the 
pipe 14. 

In certain embodiments of the present invention, a piezo-electronic pressure 

1 5 transducer may be used as one or more of the pressure sensors 15-18 and it may measure 
the unsteady (or dynamic or ac) pressure variations inside the pipe 14 by measuring the 
pressure levels inside of the pipe. In an embodiment of the present invention, the sensors 14 
comprise pressure sensors manufactured by PCB Piezotronics. In one pressure sensor there 
are integrated circuit piezoelectric voltage mode-type sensors that feature built-in 

20 microelectronic amplifiers, and convert the high-impedance charge into a low-impedance 
voltage output. Specifically, a Model 106B manufactured by PCB Piezotronics is used 
which is a high sensitivity, acceleration compensated integrated circuit piezoelectric quartz 
pressure sensor suitable for measuring low pressure acoustic phenomena in hydraulic and 
pneumatic systems. It has the unique capability to measure small pressure changes of less 

25 than 0.001 psi under high static conditions. The 106B has a 300 mV/psi sensitivity and a 
resolution of 91 dB (0.0001 psi). 

The pressure sensors incorporate a built-in MOSFET microelectronic amplifier to 
convert the high-impedance charge output into a low-impedance voltage signal. The sensor 
is powered from a constant-current somrce and can operate over long coaxial or ribbon cable 

30 without signal degradation. The low-impedance voltage signal is not affected by 

triboelectric cable noise or insulation resistance-degrading contaminants. Power to operate 
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integrated circuit piezoelectric sensors generally takes the form of a low-cost, 24 to 27 
VDC, 2 to 20 mA constant-current supply. A data acquisition system of the present 
invention may incorporate constant-current power for directly powering integrated circuit 
piezoelectric sensors. 

5 Most piezoelectric pressure sensors are constructed with either compression mode 

quartz crystals preloaded in a rigid housing, or unconstrained tourmaline crystals. These 
designs give the sensors microsecond response times and resonant frequencies in the 
himdreds of kHz, with minimal overshoot or ringing. Small diaphragm diameters ensure 
spatial resolution of narrow shock waves. 

10 The output characteristic of piezoelectric pressure sensor systems is that of an AC- 

coupled system, where repetitive signals decay until there is an equal area above and below 
the original base line. As magnitude levels of the monitored event fluctuate, the ou^ut 
remains stabilized around the base line with the positive and negative areas of the curve 
remaining equal. 

15 Furthermore the present invention contemplates that each of the pressure sensors IS- 

IS of the flow meters 10,70 may include a piezoelectric sensor 104 -107 that provides a 
piezoelectric material 1 10 to measure the unsteady pressures of the fluid/particle mixture 12 
as shown in Fig. 20. The piezoelectric material, such as the polymer, polarized 
fluoropolymer, polyvmylidene fluoride (PVDF), measures the strain induced within the 

20 process pipe 14 due to unsteady pressure variations within the process mixture 12. Strain 
within the pipe is transduced to an ou^ut voltage or current by the attached piezoelectric 
sensors 104-107. 

As best shown in Fig. 21, the PVDF material 1 10 is adhered to the outer surface of a 
steel strap 1 12 that extends around and clamps onto the outer surface of the pipe 14. The 
25 piezoelectric sensing element is typically conformal to allow complete or nearly complete 
circumferential measurement of induced strain. The sensors can be formed from PVDF 
films, co-polyraer fihns, or flexible PZT sensors, similar to that described in "Piezo Film 
Sensors technical Manual" provided by Measurement Specialties, Inc., which is 
incorporated herein by reference. The advantages of this technique are the following: 
30 1 . Non-intrusive flow rate measurements 

2. Low cost 
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3. Measurement technique requires no excitation source. Ambient flow noise is 
used as a source. 

4. Flexible piezoelectric sensors can be moxmted in a variety of configurations to 
enhance signal detection schemes. These configurations include a) co-located 

5 sensors, b) segmented sensors with opposing polarity configurations, c) wide 

sensors to enhance acoustic signal detection and minimize vortical noise 
detection, d) tailored sensor geometries to minimize sensitivity to pipe modes, e) 
differencing of sensors to eliminate acoustic noise fi'om vortical signals. 

5 . Higher Temperatures 1 40C) (co-polymers) 

10 While the present invention is capable of measuring solid particles suspended in a 

fluid, one will appreciate that other multi-phase mixtures or flows may be measured using 
an array of sensors, such as steam flow. It is further recognize the effects of dispersion on 
large solid particles in a fluid wotild be similar to large droplets of a liquid dispersed in a 
gas or air, and thus similar considerations when measuring the steam quality and droplet 

1 5 size should be addressed. 

It should be understood that any of the features, characteristics, alternatives or 
modifications described regarding a particular embodiment herein may also be applied, 
used, or incorporated with any other embodiment described herein. 

Although the invention has been described and illustrated with respect to exemplary 

20 embodiments thereof, the foregoing and various other additions and omissions may be made 
therein and thereto without departing from the spirit and scope of the present invention. 
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Claims 

What is claimed is: 

1 . An apparatus for measuring at least one parameter of a particle/fluid mixture in a 
5 pipe, comprising: 

a spatial array of at least two pressure sensors, disposed at different axial locations 
along the pipe, and each measuring an unsteady pressure within the pipe at a corresponding 
axial location, each of said sensors providing a pressure signal indicative of the unsteady 
pressure within the pipe at said axial location of a corresponding one of said sensors; and 
10 a signal processor, responsive to said pressure signals, which provides a signal 

indicative of the at least one parameter of the mixture in the pipe. 

2. The apparatus of claim 1 wherein each sensor measures an acoustic pressure and 
provides a signal indicative of an acoustic noise within the pipe. 

15 

3. The apparatus of claim I, wherein the signal processor, responsive to said pressure 
signals, provides a signal indicative of the speed of sound propagating through the mixture 
in the pipe. 

20 4. The apparatus of claim 3 wherein said signal processor comprises logic which 
calculates a speed at which sound propagates along said spatial array. 

5. The apparatus of claim 3 wherein said signal processor comprises logic which 
calculates a frequency based signal for each of said acoustic pressure signals. 

25 

6- The apparatus of claim 4 wherein said acoustic pressure signals each comprise a 
frequency based signal and wherein said signal processor comprises logic which calculates 
a ratio of two of said frequency based signals. 

30 7. The apparatus of claim 1 comprising at least three of said sensors. 
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8. The apparatus of claim 3 wherein the signal processor comprises logic which 
calculates a fluid composition of the mixture in the pipe. 

9. - The apparatus of claim 1 wherein at least one of said pressure sensors measures a 
5 circmnference-averaged pressure at said axial location of said senson 

10. The apparatus of claim 9 wherein at least one of said pressure sensors includes a 
piezoelectric sheet material. 

10 11. The apparatus of claim 1 wherein the piezoelectric sheet material is polarized 
fluoropolymer, polyvinylidene fluoride (PVDF). 

12. The apparatus of claim 1 wherein at least one of said pressure sensors measures 
strain on the pipe. 

15 

13. The apparatus of claim 3 wherein the frequency based sound speed is determined 
utilizing a dispersion model to determine the at least one parameter of the mixtm-e. 

14. The apparatus of claim 3 wherein the array of acoustic sensors are spaced 

20 sufficiently such that the entire length of the array is at least a significant fraction of the 
measured wavelength of the acoustic waves being measured. 

15. The apparatus of claim 1 wherein at least one sensor provides a first filter which 
measures a vortical pressure field at a first axial location along the pipe and provides a first 

25 pressure signal indicative of said vortical pressure field; and 

at least a second sensor provides a second filter which measures said vortical 
pressure field at a second axial location along the pipe and provides a second pressure 
signal indicative of said vortical pressure field. 
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16. The apparatus of claim 15, wherein the signal processor, responsive to said first and 
second pressure signals, provides a velocity signal indicative of a velocity of the said 
vortical pressure field moving in the pipe. 

17. The apparati^s of claim 16 wherein said velocity signal is related to a velocity of said 
fluid moving in said pipe. 

1 8. The apparatus of claim 1 5, wherein said first and said second filters filter out 
wavelengths associated with an acoustic pressure field and passes wavelengths associated 
with said vortical pressure field. 

19. The apparatus of claim 18, wherein said first filter comprises a first spatial filter; and 
said second filter comprises a second spatial filter. 

20. The apparatus of claim 19, wherein: 

said first spatial filter comprises at least a first and a second unsteady pressure 
sensors di^osed a predetermined first distance apart from each other; and 

said second spatial filter comprises at least a third and a fourth unsteady pressure 
sensors disposed a predetermined second distance apart from each other. 

21. The apparatus of claim IS wherein said signal processor comprises logic which 
calculates a cross-correlation between said first and said second inhomogeneoxis pressure 
signals and provides a time delay signal indicative of the time it takes for said vortical 
pressure field to move firom said first location to said second location. 

22. The apparatus of claim 21 wherein said signal processor comprises logic responsive 
to said time delay signal which provides an inhomogeneous velocity signal indicative of the 
velocity of said vortical pressure field moving in said pipe. 
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23. The apparatus of claim 21 wherein said signal processor comprises logic responsive 
to said time delay signal which provides said velocity signal indicative of the velocity of 
said fluid moving in said pipe. 

5 24. A method for measuring at least one parameter of a particle/fluid mixture in a pipe, 
said method comprising: 

measuring unsteady pressures within the pipe at at least two predetermined axial 
measurement locations along the pipe to provide a pressure signal indicative of the unsteady 
pressure within the pipe at each of the at least two predetermined axial measurement 

10 locations; and 

calculating the at least one parameter of the particle/fluid mixture in the pipe using 
the tmsteady pressure measured at the axial measurement locations. 

25. The method of claim 24 wherein the measured unsteady pressures are acoustic 
15 pressures to provide a signal indicative of an acoustic noise within the pipe. 

26. The method of claim 25, wherein the calculating the at least one parameter uses an 
acoustic pressure to calculate a speed of sound propagating in the pipe, 

20 27. The method of claim 24 wherein measuring iinsteady pressures comprises: 

measuring a vortical pressure field at a first location along the pipe and providing a 
first vortical pressure signal indicative of said vortical pressure field; and 

measuring said vortical pressure field at a second location along the pipe and 
providing a second vortical pressure signal indicative of said vortical pressure field, said 
25 first and said second locations heing an axial distance apart 

28. The method of claim 27, wherein the calculating the at least one parameter uses the 
first and second vortical pressure signals to calculate the velocity of the mixture in the pipe. 

30 
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